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ABSTRACT 


A mathematical model of a vehicle being brought to rest in 
such a manner as to cause a skid is devised and simulated on an 
analog computer. This model includes the effects of the brake 
line pressure upon the braking of the vehicle and also the effects 
of the coefficients of friction between the tire and road surface 
and between the brake shoe and brake drum. Several controls to 
correct for the occurrence of a skid are simulated on the analog 
computer and a detailed analysis of the braking system under the 
influence of each control is performed. Lastly, a scheme for 
implementing the control which minimizes the stop time and the 


degree of skid is presented. 
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1. INTRODUCTION 


Skidding of a vehicle due to an excessive application of 
brake force for the road conditions prevailing causes many 
thousands of accidents every year. Rapid and heavy application 
of brakes in a panic stop situation often causes wheel lock up 
and subsequent loss of control of the vehicle. This also occurs 
in wet or icy road conditions. 

The minimum stopping distance possible for a vehicle is 
determined by the coefficient of friction between the tire and 
the road surface. When a vehicle is just on the verge of a skid 
the frictional force between the tire and road surface is at its 
greatest value, and consequently the deceleration of the vehicle 
is at a maximum, i.e., equal to M8: where ae is the coefficient 
of static friction between the tire and the road surface and g 
is the deceleration of gravity (see Appendix A for a proof of 

(7: 
this statement). When skidding occurs the tire slides along the 
road surface. This causes a reduction in the coefficient of 
friction and a corresponding reduction in the frictional force 
between the tire and road surface, thus decreasing the magnitude 
of the vehicular deceleration. The optimum braking situation 
would thus be attained if the vehicle could be kept just at the 
verge of a skid. If a skid could be detected exactly at the 
instant of occurrence and a compensating force were automatically 
applied to eliminate the skid before wheel lock up occurred, then 
nearly ideal braking conditions could be realized under all road 


conditions. 


Attempts have been made previously to develop "anti-skid" 
devices, and some such devices are presently in use on certain 
aircraft and many trucks. Previous efforts in this direction 
include the Dunlop System of England, the Maxaret Unit, the Anti 
Skid Device, the Westinghouse Decelostat Controller, the Hydro- 
Aire Hytrol system, and the Anti-Skid Braking system of Russia's 
Tu - 104 Jet Liner.[1] These devices all work satisfactorily to 
a greater or lesser extent, but most have the drawback that 
operation depends upon the sensing of either a complete wheel 
lock up or the attainment of a pre-determined rate of wheel 
deceleration, at which time a bang-bang type of control is initi- 
ated. Additionally, many of these devices require expensive, 
heavy installations, rendering them largely unsuitable for eco- 
nomical installation and operation on commercial vehicles. 

This thesis presents an analysis and evaluation of several 
electrically actuated alternate proposals designed to automatically 
compensate for an overbraking situation in order to minimize the 
degree of skidding and more nearly optimize the braking of a vehicle. 
A mathematical model for the system was developed and simulated on 
an analog computer, and three types of skid correcting controls 
were devised and tested. One form of control which linearly 
decreases the net force applied to the brake shoe until the skid 
is eliminated and then permits a linear er easc of the net braking 
force until a skid once again occurs was tested. This is called 
the double ramp control. Two forms of bang-bang controls were also 
tested. Finally a triple ramp control which continually "samples" 
the road surface to determine the optimum brake force was devised 


and tested. 
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Complete data were obtained for each type of control and a 
detailed analysis was made of the behavior of the braking system 
under the influence of each control. 

Lastly, a scheme for implementing the triple ramp control 
(which proved to be superior to any others devised) is advanced, 


although the physical system was not actually built. 


dh 


2. DESCRIPTION OF THE PROBLEM 


For the purposes of this paper, a skid will be defined as 
that condition which exists when the actual linear velocity of 
the frame of the vehicle in question is greater than the product 
of the angular velocity of the wheel and the radius of the wheel. 


In terms of the symbology used throughout this paper 


e > 0 when vrRyo (2) 
and e = 0 when v = RIO; (3) 


where v represents the linear velocity of the frame, R, represents 


L 
the radius of the wheel and is the angular velocity of the wheel. 
The term e corresponds to an error signal, the magnitude of which 

is a measure of the degree of skidding present. The product Rw 
will be referred to as the velocity measured by the rotation of the 
wheel. When e = O there is no skidding and when e > O the vehicle 

is moving faster than a point on the rim of the wheel, i.e., the 
wheel is tending to lock up and a skid is occurring. The condition 
when Rw V would cccur when the wheels are spinning, as, for example, 
when trying to accelerate on an icy road. 

A proposed means of sensing the linear velocity of the vehicular 
frame is to integrate the output from a linear accelerometer attached 
to the axle housing of the vehicle (to minimize effects of vehicle 
tilt during braking), and to compare the integrator output with the 


velocity measured by the rotation of the wheel. The initial condition 
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input to the integrator would be the velocity measured by the 
rotation of the wheel at the instant of application of the brake 
force. 

When an error signal occurs, indicating skidding is present, 
an automatic corrective force (FD would be generated to reduce 
the total force F applied to the brakes. Fin represents the force 
applied by the vehicle operator to the brakes and is considered 
constant for all subsequent analyses. 

These concepts are illustrated in Fig. 1 where a block diagram 


of the proposed system is eh Gahae 


Vo= R,Wo 


east 


LINEAR nag | ae 


Vig V [oom 
ACcELER- | = f ji —->@-5 = fF, 
OMETER ie 


= INTEGRATOR | 


Fig. 1. Block diagram of the proposed braking system. 


From elementary physics it is known that the coefficient of 
friction between two surfaces is dependent upon the relative velocity 
between the surfaces. A sketch of the behavior of the coefficient 
of friction as a function of the relative velocity between the two 


; nee [8] 
surfaces is shown in Fig. 2. 
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Fig. 2. Coefficient of friction as a function of relative velocity. 


Just before relative motion begins, the coefficient of friction 
between two surfaces is at its maximum value. Prior to reaching 
this maximum value, the coefficient of friction is a variable that 
adjusts itself in such a manner so that the force tending to cause 
a relative motion is just compensated by the reactive force of 
friction (this is before any relative motion takes place). After 
sliding motion starts (as between a brake shoe and drum, or between 
a tire and road surface for a skidding situation) the coefficient of 
friction # varies as shown. 


A piecewise linear approximation to this behavior is shown in 


Page 3 4 





ee 
Ve RELATIVE 
VELOCITY 


Fig. 3. Piecewise linear approximation to the behavior 
Of thé coeGfEI client Of friction. 
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The term be corresponds to the coefficient of static friction 
which is a maximum just before relative motion begins, and which is 
always greater than the coefficient of kinetic or sliding friction. 


The equation for this piecewise linear approximation is 


fh = Keo , for Oa Ya Ve (4) 


Mee Me —K, (Vee) , For Ve € VE(Ve + ) (5) 


The above relationship will be used to represent the coefficient 
of friction between the brake drum and brake shoe of the moving vehicle. 
The break point, ee for the coefficient of friction between the tire 
and the road surface will be assumed to be high enough so that for 
normal speeds the coefficient of kinetic friction is some constant 
ae As will be seen, this is not a restrictive assumption, since the 
purpose of the corrective forces to be applied is to minimize the degree 
of skid experienced. Thus the relative velocity between the tire and 
road surface would be small enough to warrant considering the coefficient 
of kinetic friction to be a constant. 


Consider the idealized brake drum and shoe shown in Fig. 4. 





C) 
> BRAKE 
TIRE ee ie ce 
> ‘ Sua x , eae ein — es 


KRa 
! Se 


{ SS 
rR, : e/) | Ty, =Net Porque 
ayes “A / R,= Radius of Wheel (Tire) 
. Re= Radius of Brake Drum 
= ear 


Fe 


Fig. 4. Diagram of brake drum and shoe. 
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In Fig. 4 the wheel is assumed to be rotating in a clock-~ 
wise direction with angular velocity W and angular acceleration 
W with the positive sense as shown. F represents the total (net) 
force acting on the brake shoe, and is the difference between Fn 
(applied by the brake pedal) and oe (the corrective force). Fy 
is the force of friction, the value of which depends upon whether 


a skid is present. The equations describing this system are pre- 


sented below. 


i aren (6) 
Net torque = T =R, Fy - R,F, (7) 
Fe = fags” yutTor Ee. =o (8) 

fermg ») for e>O (9) 


where m = the weight of the vehicle acting through the axis of 


the wheel divided by the acceleration of gravity g (= 16 fejecc a 


> 


[A= feo~ K/ Re (W- We), for we We (10) 


fro ) a ea CC) Coe. Gl) 


i.e., | is the coefficient of friction between the brake drum and 


brake shoe (Fig. 3). Also, is the coefficient of friction between 


BK 


the tire and the road when a skid is present. 


Let ie Re= K, , so that Eq. (10) becomes 
[= e- K,(W-We) for WP We (12) 
[0 } for Of We (13) 


and Eq. (7) becomes 
T=ReFe Rk, Fe = —I a (14) 


where J = the moment of inertia of the wheel. 
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Consider the case where e = 0, i.e., skidding is absent, so 


Fee —-mv & Meg , (15) 
with the equality holding when the system is just at the verge of 
a skid. The equation of motion, (14), becomes 

-Ta = R,mv+ReFpu . (16) 
Since e = QO, then 

eco (17) 
Substituting (17) into (16) yields 

~Jo = Ri mo +RoFp . (18) 


Solving (18) for ™ yields 


ee hae: (19) 
J+ R2'm 


When e @ 0, i.e., skidding is present, and 


eal kk mg , where [er fs | (20) 
then (14) becomes 

=Jiw=RePA- R,Anmg ; a 
SO 

a Rn mg _ —— ar (22) 


Note that for each of the relations for W the value of uw is 
a function of W, In the case where W ee then » = ee and the 


equations above apply with » replaced by Mee When W 2 Dees then 


Bafeo — K,Cw- We). _) 


The complete set of equations is indicated below: 


When e = 0, W< We 3 


O= -~Re F 2 (24) 
T+ R= rn 
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When e > QO, Cea CO 


C) = RK, fA kk cee _ Re fe, a (25) 


When e = 0, © 2We: 


w= ReF K, g —~ RoFKa , (26) 
J+Rym JT+Rem 


When e 7? 0, O> Wc: 


= ReaEK, ~ — ReF ke + Rifxm@ , (27) 
FT J : 


In the last two relations, 
Ke Mot Ki we . 

The analog computer simulation of this system of equations is 
shown in Fig. 5. Note that Appendix B contains a brief description 
of the analog computer symbols used in this paper and an explanation 
of the logic elements used to control various portions of the simu- 
ieee 

To determine the equations governing the forward motion of the 


vehicle, recall that the maximum possible deceleration for a given 


coefficient of friction w' between the tire and road surface is 


ma ee a } (28) 
where f= Ms of KK. 


Thus, when e = 0, i.e., no skidding, 
-x=-Riwof Msg ; (29) 
where be = the coefficient of static friction between the tire and 


the road surface. Also, before a skid occurs, 


x =R w= y, (30) 
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After a skid occurs, 


x = eG J (31) 
where be = the coefficient of kinetic friction between the tire 


and road surface which is considered constant for this simulation, 


ae ee 
Thus, when e > 0, 
Ne gt + Vo 3 
where Ve is the velocity of the vehicle when the skid begins. 
The analog computer implementation of these equations and 
ideas is presented in Fig. 6, which provides a simulation for V 


under conditions of skidding or skidless braking. 








This amplifier begins integrating 
when €>0 or R,w <-Msg- 
aa Otherwise it i$ ihn reset and 


the initial Conditiong appear 
at the output, 


Fig. 6. Analog computer simulation of vehicular velocity. 


The composite analog computer set up for both the brake drum 
and shoe and the vehicular velocity is shown in Fig. 7. This is 
intended to represent a conventional braking system in existence 


on present day vehicles and no corrective force is included in 


this diagram. 
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3. PROPOSED SKID CORRECTING CONTROLS 


Any form of corrective control must be initiated by the 
occurrence of a skid, which in turn is indicated by the presence 
of a positive error signal. The purpose of the corrective control 
must be to lessen the net force F applied to the brake shoe when 
a Skid is sensed. When the skid has ceased, the corrective force 
F. must decrease to permit a reapplication of the full braking 
force to the brake shoe. Thus the value of F. will depend upon 
whether a skid is present or not. 

The first form of control to be tried is called a double ramp 


control, and a typical shape for the control is shown in Fig. 8. 





C70 


Fig. 8. Double ramp control. 


An examination of Fig. 8 indicates that the maximum value 
attained by B is Bea (the force applied by the vehicle operator). 
Since F = con - a , and since there can never physically be a 
negative force on the brake shoe, BP. = Faw The value of K, is 
not necessarily equal to Kye 


An analog computer simulation of the double ramp control is 


shown in Pie. 9. 
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Fig. 9. Analog computer simulation of the double ramp control. 


The integrator indicated in Fig. 9 may be in the compute mode 
for the entire problem simulation, due to the presence of the feed- 
back limiter at the output of the integrator which prevents PF. 
from ever going negative. 

A bang-bang control was tried for the second form of corrective 
force. Since any physically realizable bang-bang control would have 
either a dead zone or a time delay, one control with a dead zone 
and one with a time delay were tried. The expression "dead zone" 
is meant to imply that there is some positive value of error volt- 
age above which the corrective force is full on until the error is 
once again reduced to zero. The analog computer implementation of 
the bang-bang control with a dead zone is shown in Fig. 10, where @ 


is the value of the dead zone voltage. 





Fig. 10. Analog computer implementation of the 
bang-bang control with a dead zone. 
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The expression "time delay'' as used here implies that the 
error exists for some time Ty before corrective action is initiated. 
In order to simulate the time delay, a Wood fourth-order approximant [2] 
time lag circuit was set up to delay the error voltage by one second 
(real time delay), which in problem time corresponded to a delay time 
of 0.1 second. The analog computer diagram for the bang-bang control 


with a real time delay, with the potentiometer settings needed to 


produce a one second real time delay, is shown in Fig. ll. 





On Te 
nw Ut-) > 0 
ft 


ae —_ 


Fig. 11. Analog computer implementation of the 
bang-bang control with a real time delay. 


The third form of control to be simulated is called the triple 
ramp control, and a typical plot of corrective force F. versus time 
plus the complementary plot of F (= Fen = FO) versus time is shown 


in Figw 12. 
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Fig. 12. Triple ramp control. 


The triple ramp control is of the following form. F. = 0 
until a skid occurs. When a skid occurs, e > O and FO begins 


rising rapidly with a slope K, to a maximum value of Fn or until 


5 
e = 0. When e = QO, F. then begins decreasing with a slope equal 

to - Ky until e > 0 once more. Up to this point the control is 
identical to the double ramp control. However, the second time 

that e > 0 a separate monitoring device with a "memory" (a Mylar 
capacitor, for example) holds the voltage at the value corresponding 


to the corrective force at which the skid occurred, 


At the instant when e > 0 for the second time in Fig. 12, F 


c 
would once again increase rapidly with slope K, until e = 0. Then 
Ee would again begin decreasing with slope - K,, until it is just 


above the corrective force at which skidding occurred previously. 
Then Be would be allowed to decrease more slowly with a slope - K. 
(thus permitting a gradual increase in F until a skid once again 


occurs). The cycle then repeats itself. 
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The design of this particular type of control was motivated 
by the considerations expressed in the Introduction, i.e., that 
the optimum braking situation would be realized if the vehicle could 
be kept just at the verge of a skid. The storage of the value of 
F. at which a skid occurs endows the control with a form of "intel- 
ligence' which permits the net force applied to the brake drum to 
be just less than the optimum force. The gradual increase in F 
permitted by the third portion of the control (with a slope of - K.) 
is necessary to permit the control to adapt to improving (i.e., less 
slippery) road conditions. The first cycle of the control, where 
ES increases rapidly with a slope of K, and then decreases with a 
slope of - Ky until the skid peetee once again, is necessary to 
permit the control to "test'' the road surface for the net force 
required to just cause a skid. The analog computer implementation 
of this control is shown in Fig. 13. 

Integrator No. 1 provides the ramp with a slope of K,, - Kyo 
or - K. as appropriate. This integrator is in the compute mode for 
the entire simulation time and F. is kept from going negative by the 
lower bound of the feedback limiter. The upper bound on the limiter, 
and conseguently on Fo is determined by Peat Potentiometer No. l, 
labeled C.R./10, determines the value of Fat which the slope 
changes from - Ky, to - Ko3 the? rabnosok Fo) Feta) (see Fig. 12) 
is termed the Correction Ratio (C.R.). It is this ratio which 
determines the setting of potentiometer No. 1. When integrator No. 2 
is in the track (reset) mode switch No. 1 is closed (FF = 1) and we 
desire the slope of F. to be either K, or - Ky, (but not ~ K); since 


this is the time period between the first and second occurrences of 


the positive error signal. Thus to positively prevent relay No. l 
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Fig. 13. Analog computer implementation of the triple 


ramp control. 


Lap 


from throwing during this phase of the cycle potentiometer No. 2 
is set lower than potentiometer No. 3 to insure that C = 1 and 
C = 0, where 1 and O correspond to a logical true and false, 
respectively. 
The flip flop (the output of which is denoted by FF and FF) 
is an R - S flip flop, and is used to control the mode of integrator 


No. 2 as well as the state of switches No. 1 and 2, 
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4. THE UNCONTROLLED SITUATION 


The analog computer simulation of the uncorrected braking 
of the vehicle (shown in Fig. 7) was set up on the COMCOR Ci 


5000 Analog Computer, using the following arbitrary parameters: 


Ky = 0.01 
. 60 rad./sec. 
Lo = 0.8 
Thus wu = 0.8 - 0.01 (Ww - 60), we = 
and Ww = 0.8, W< Wo» 
where WU = coefficient of friction between the brake drum 
and shoe. 
bo = 0.75 = static friction between tire and road surface. 
by = 0.5 ==> kinetic friction between tire and road surface. 
—_ 75 £t./sec. = initial vehicular velocity. 
a= 75 rad/sec, = initial rotational velocity of the wheel. 
Ry ee aia 
Ro =O Ee. 


J =225 ft. Ib. weer 


M = 93.75 slugs 


(Weight = mg = 2000 lbs. = weight of vehicle acting 
through axis of wheel). 


A 7000 1b. input force (corresponding to FP. was simulated as 
a 70 volt step input occurring at time = O sec. 
Six different variables were plotted for each run on a Brush 
Instruments Mark 200 (6 channel) Recorder. They were 
W ==> the angular acceleration of the wheel, in Pad seca 


W ==> the angular velocity of the wheel, in rad/sec. 


Zo 


V = the actual linear velocity of the vehicle, 
int fey seer 


F = > the net force applied to the brake shoe 
(= Pn - BO» in lbse 


e => the error signal (=v - Rj), inetp/sec. 
B =? the corrective force, in lbs. 

The uncorrected system response of the vehicle, i.e., the 
response obtained when no corrective force is applied, is shown in 
Fig. 14. From examination of the uncorrected system plot it is 
seen that the wheel locks in 0.96 sec. and skidding occurs 0.16 
sec. after application of the braking force. The skid is present 
until the vehicle comes to rest 4.64 sec. after the force is initially 
applied. 

Measuring the slope of the velocity output plot gives an 
average deceleration of 16 Felace: (excluding the first 0.1 sec. 
of the problem, where skidding does not occur). This is expected, 


since with ., = 0.5 the deceleration during a skid is 16 Eee 


K 
(=9025-x 32 rt cecmane 

To completely understand the shapes of the various plots 
displayed in Fig. 14, recall that the coefficient of friction between 
the brake drum and brake shoe is of the general shape indicated in 
Fig. 3, with a break point at w= 60 rad/sec. The initial rotational 
velocity of the wheel was a= 75 rad/sec., and thus at the time of 
simulated application of the brakes the coefficient of friction iu 
between the brake drum is 

u = 0.8 - O.CIM/S —360)5— 0765. 
The value of wu will increase as the relative velocity between 


the brake shoe and brake drum decreases, until it has attained a 


maximum value of 0.8. After the problem begins there is a time lapse 
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Fig, 14, Uncorrected braking situation. 
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of 0.16 sec. before skidding starts. During this initial 0.16 
Sec. is increasing since the vehicle is slowing down, However, 
uw is sufficiently small during this time interval so that the net 


torque T (= Ro Fu - RLF ) acting on the wheel is not sufficient to 


£ 
cause a skid. When the velocity has decreased to 70 miles/hr., 


and w is down to 70 rad/sec. (since a skid has not occurred yet), 


the value for » has increased to 


» = 0.8 - 0.01 (70 - 60) = 0.7. 

This value for & is sufficiently large to cause the net torque 
to be just equal to that value necessary to cause a skid. At the 
instant a skid occurs the coefficient of friction between the tire 
and road surface drops abruptly from b= O57 DEO be = OR" Fe also 
decreases abruptly due to the occurrence of the skid and the net 
torque increases rapidly, causing the wheel to decelerate rapidly 
and lock after 0.96 sec. total elapsed time. 

Note that the plot for the radial acceleration of the wheel W 
indicates that at time = 07 the acceleration drops sharply to a 
negative value (- 23 Paevicce and then slowly becomes more negative 
until the time when a skid occurs. During this time the wheel is 
still slowing down gradually, and the deceleration (or negative 
acceleration) is increasing slowly as the value of » (between the 
brake shoe and brake drum) increases. At the point when the value 
of W becomes just less than ~ 24 ral eec the skid occurs, and W) 
abruptly becomes much more negative (- 76 Pe ceene) and decreases 
rapidly until it reaches a constant - 94 ped eecne The significant 
factor to be noted here is that the skid occurs just when W = 


- 24 rad/sec.’. This corresponds to a value of V = - 24 fe) eee 


(since Ry = 1 ft., and a skid has not occurred yet). This is in 


SZ 


agreement with the previous statement that the maximum deceleration 
possible without a skid occurring is 
oa Oey Dek Oe ft/sec.” = - 24 Pe een 

The leveling off of the w plot at - 94 Pad eece occurs at 
the time when the rotational velocity “ crosses the 60 rad/sec. 
mark which is the value established for the critical value of 
rotational velocity W below which BP becomes a constant = 0.8. 
Since & is constant and 

Fe = Uemg 

then the net torque T is 


IT = Ry Fu - Rib ng = a constant, and W is thus a 
constant. 


The magnitude of the error rises to a maximum value and then 
decreases to zero when the vehicle has completely stopped. This 
behavior is expected when it is noted that the error is merely the 
difference between the linear velocity and the rotational velocity 
of the wheel. 

One other comment concerning the response of the uncorrected 
System is appropriate. Inspection of the vehicular velocity v plot 
reveals that a steeper slope is present on the plot prior to the 
beginning of a skid. After the skid begins the slope becomes less 
steep and remains constant until the vehicle comes to rest. The 
deceleration of the vehicle is thus greater before the skid begins 
than it is after skidding occurs, which is in accordance with the 
expected results. 

The important fact to notice is that during a large part of 


the path the vehicle is in what would be an uncontrolled skid, as 
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shown by the large magnitude of the error. This uncontrolled skid 

is a matter of as much concern as the lengthened time to stop the 
vehicle. This lengthened time is caused by the effect of the 

smaller coefficient of friction which acts between the tire and 

the road surface when skidding is occurring. The optimum time of 
stopping, assuming braking takes place just at the verge of skidding, 
1s 


75 ft/sec. 


24 ft/sec® Bo co 


t optimum = 
jae ‘ ; : 
where 24 ft/sec. is the maximum possible deceleration uSing the 


values assumed Eor EHS Simulation. 
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5. THE DOUBLE RAMP CONTROL 


The various proposed controls described in Section 3 were 
implemented in an attempt to obtain a more nearly optimal braking 
Situation. The first of these corrective forces was the double 
ramp control shown in Fig. 8. A total of 21 separate runs were 


made with this form of control. K, was initially held constant 


A 
at 40000 lbs/sec., and Ky was varied from 5000 lbs/sec. up to 
60000 lbs/sec. Then K. was held constant at 60000 lbs/sec. and 


3 


Ky varied from 5000 lbs/sec. to 50000 lbs/sec. Finally, two 


runs were made with K, and K, set at 80000 lbs/sec. and at 100000 
lbs/sec. Representative samples of the behavior to be expected 
from this form of control are shown in Fig. 15 where KR, = 60000 


lbs/sec. and Kaas 35000 lbs/sec., and in Fig. 16 where K, = 


100000 lbs/sec. and K, = 100000 lbs/sec. 


2D 
An examination of Fig. 15 indicates that for the specific values 
of K, and K, chosen (i.e., 60000 lbs/sec. and 35000 lbs/sec. respec- 
tively) the stop time has been increased to 5.47 sec., whereas the 
stop time under controlled conditions was only 4.64 sec. However, 
the total time during which a skid is occurring has been reduced to 
2.08 sec. and the magnitude of the skid has been reduced to approx- 
imately 2 ft/sec. during the entire braking cycle. These figures 
are contrasted to a total skid time of 4.47 sec. and a maximum magni- 
tude of skid of 68 ft/sec. for the uncorrected situation. 
When the brakes are initially applied at time = 0 sec. there is 
a time lapse of 0.16 sec. before a skid occurs. This is exactly the 


same situation that occurs in the uncorrected braking and is accounted 


for by the same reasoning. After 0.16 sec., a skid commences. This 
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Fig. 16. Double ramp control response, K., = 100000 lbs/sec., 


Ky, = 100000 lbs/sec. 


is indicated in Fig. 15 by the positive going error signal and also 
by the abrupt drop of the W plot from - 24 Pe eer to approximately 
- 78 pedi seman Recall that when skidding occurs the net torque 
acting to slow the wheel is suddenly increased, thus the wheel 

starts to slow down rapidly, as seen from inspection of the Ww 

versus time plot at the time of the initial occurrence of the first 
skid. However, the corrective force immediately begins acting to 
decrease the net force F (= Pn ~ PF.) acting on the brake shoe, 

which consequently decreases the net torque T (= RAFU - Ri bmg) » 

and the wheel begins to speed up to "catch up" with the velocity of 
the vehicle. This accounts for the sudden positive slope of the Ww 
plot almost immediately after an error signal occurs, and also for 
the upswing on the W versus time plot after the initial drop of w 
with each occurrence of a skid. The corrective force increases until 
it equals Fen and/or e = 0 once more, at which time it begins decreasing 
towards zero until a skid occurs again. 

An expanded plot of W and v versus time for two skid cycles is 
shown in Fig. 17. To obtain this plot the equations of motion of the 
system were solved numerically by slide rule with the value of K, = 
5000 lbs/sec. and Ky, = 8000 ibafeee oe the simplifying assumption 
was made that the coefficient of friction between the brake drum and 
brake shoe was constant at 0.8. While these values for K, and K, 
and the simplifying assumption for @ were not used in the analog 
computer runs, the shape of the curves is similar to the shape of 
those obtained on the analog computer simulation. A brief explanation 


of Fig. 17 will facilitate understanding the curves shown in Fig. 15 


and Fig. 16. 
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Referring to Fig. 17, note that a skid begins. at time = O° 
sec. This is due to the simplifying assumption of a constant 
UL = 0.8 over the entire velocity range of interest, and in no way 
affects the qualitative aspects of the curve shown. The value of 


W (or R.W as plotted in Fig. 17) decreases initially and then rises 


it 
to intersect the velocity V plot (which is decreasing linearly as 
the vehicle skids). For time less than the point marked Gr the 
slope of the W curve would be positive, as can be seen by mental 
differentiation of the Rw curve. This is in agreement with the w) 
curve of Fig. 15 and Fig. 16. At the point marked Gis however, 
the Rw curve now begins increasing linearly, and the W curve would 
have a constant value, i.e., would level off at some positive value 
corresponding to the slope of the segment from (1) to Oe This 
leveling off of dune plot is observed also in the analog computer 
results of Fig. 15 and Fig. 16. The time period when the leveling 
off occurs corresponds to the region between () and (2) of the 
inset plot of Bao 1.e., where PF. has attained its maximum value 

and the net torque is now a constant. Since the net torque is 
constant, the value of W is also a constant, until the skid is 
completely eliminated. 

At point (2) the skid is eliminated and F. begins dropping to 
permit reapplication of the full brake force Fea’ However, the net 
force F increases gradually from zero with a slope K,, and for a 
Significant time after the elimination of the skid the velocity 
is decreasing at an unacceptably slow rate. This is seen by an 


inspection of the region of the curve of Fig. 1/7 from point (2) to 


@- At point 6) the deceleration has finally increased to that 
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value of deceleration experienced in a skid, i.e., - 16 OES. 
and at point (4) it is sufficiently great to cause a skid to occur 
once more. The process then repeats itself. 

The rather detailed qualitative analysis undertaken above serves 
to point up the major defects inherent in the double ramp form of 


control. If the parameters K 


3 and Ky are Chosen so as to maximize 


the time when a skid is not occurring, e.g., to maximize the time 
from point (2) to point (4) of Fig. 17, then the time from point (2) 
to (3) will also be maximized. From (2) to ©) the velocity is 
decreasing at less than the optimum rate, and even at less than 
the deceleration experienced when in a skid. Thus it would seem 
that by minimizing the skid time with the double ramp control, we 
are inadvertently maximizing the actual time to stop. This conclusion 
is borne out by the results of Fig. 18 and Fig. 19, where the skid 
times and stop times for the various values of K, and Ky, simulated 
on the analog computer are plotted. 

Figure 20 is a plot of the magnitude of the skid obServed versus 
K,> with K,, fixed at 40000 lb/sec. The magnitude of the skid after 
it has obtained a constant value is the value that is plotted here. 
It is seen that the magnitude of the skid is very dependent upon the 
value of K,; which is an expected result since the larger the value 
of Ky» the more rapidly is the skid corrected. This also accounts 


for the decrease in the skid time with increasing K, observed in Fig. 


3 
18. Experimental results indicated that the magnitude of the skid is 
unaffected by the value of K,, which is natural since the downward 


going portion of the double ramp control is only acting when no skid 


is present. The stop time is increased, however, with increasing Ky 
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Fig. 18. Stop time and skid time versus K, for the 


double ramp control. 
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Fig. 20. Magnitude of skid versus K. for the double 
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ramp control, 
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or Ky Since, as mentioned previously, the greater the proportion 
of the total stop time that is free from a skid, then the greater 
will be the actual time when the deceleration is less than that 
experienced when a skid is present. 

Even though the double ramp control has serious defects, it 
nevertheless represents a decided improvement over the completely 
uncontrolled braking represented in Fig. 14. Since under many road 
conditions it is the loss of control which occurs when a skid is 
present that is actually more dangerous than the increase in the 
stop time, then the double ramp control would have certain merits. 
The double ramp control is capable of minimizing the skid time and 
the magnitude of the skid, at a cost of slightly increased stopping 
distance. 

If it is assumed that the stop time, the skid time, and the 
degree of skid experienced are of equal importance in judging the 
overall goodness of a particular type of control, then a performance 
factor (P.F.) can be defined for the uncorrected case, the double 
ramp control, and for subsequent controls to be tried. This per- 
formance factor is defined as the reciprocal of the product of the 
actual stop time minus the optimum stop time, the skid time, and 
the degree of skid experienced. The degree of skid will be measured 
after the velocity has decreased below 60 ft/sec., since that is the 
velocity below which the value for uw (between the brake drum and 
brake shoe) is a constant. The optimum stop time is 3.13 sec. 

This definition of a performance factor must not be taken to 
be the absolute criterion for judging a control and can only be 


used as a rough guide to determine the relative merits of different 
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types of controls. It cannot, for example, be used to determine 
the proper setting of control parameters to give the ideal stopping 
conditions, since the optimum settings for one set of road conditions 
might very possibly be different from the optimum settings for another 
set of road conditions. Also, no account is taken, in the definition 
of a performance factor, of the fact that a slight degree of skidding 
might be much less significant than a larger degree of skidding, or 
that an increase in the stop time might be permissible under some 
circumstances but entirely unacceptable under others. Nevertheless, 
the definition given above is convenient and useful, within limits. 

For the uncorrected situation of Fig. 14, a performance factor 
can be computed in order to provide a basis for comparison. The 
integral of the error curve, rather than the product of the error 
times the skid time, is used, multiplied by the total stop time minus 
the optimum stop time. This gives a performance factor of 0.00214 
for the completely uncontrolled situation. 

For the double ramp control, the values of the performance 
factor range from 0.0154 to 0.298, and the variation of P.F. with 


a change in K, is shown in Fig. 21. This figure is presented only 


3 


to indicate the trend of P.F. for the double ramp control and not as 


an aid to selecting the "best" value for K However, based on the 


a4 
indication given by the performance factors for the double ramp 


control, the double ramp control provides a definite overall perform- 


ance of the system, when compared to the uncontrolled situation. 
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Fig. 21. Performance factor versus K., for the double 


ramp control. 
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6. THE BANG-BANG CONTROLS. 


The bang-bang control with a dead zone was simulated on the 
analog computer. The dead zone was arbitrarily set at 7.5 ft/sec., 
which is rather high for a physical system, but which nevertheless 
serves to illustrate the behavior of the system under this form 
of control. The response obtained is shown in Fig. 22. 

The stop time is seen to be 4.50 sec., which is slightly less 
than that obtained for the completely uncontrolled situation. This 
is a predictable result, since for almost the entire cycle the 
vehicle is in a "controlled" skid, where the term "controlled" 
implies that the value of the skid is never permitted to exceed a 
certain value. The reason for the slightly lower stop time than 
that observed for the uncontrolled situation is that for a small, 
non-zero time interval in the region where F. drops to zero and F 
rises to Fn? the value of the skid is essentially zero (or at any 
rate less than the threshold established for detection of an error 
on the basic system of Fig. 7). While the skid is this small, the 
coefficient of friction controlling the deceleration of the vehicle 
is be (which is greater than bie)» and the deceleration is thus 
greater for this small time interval. Since this phenomenon occurs 
every time the skid is corrected, the cumulative effect over the 
entire stopping distance is to decrease the stop time slightly. 

The major disadvantage to this form of control is the fact 
that a skid, even though of a small value, is almost always present 
until the vehicle comes to a complete stop. This also implies that- 
the shortest possible stop time is nearly 4.64 sec., i.e., the same 


time required to stop when no corrective force is present. However, 
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the fact that the skid is controlled, and can theoretically be 
reduced to an arbitrarily small positive value (by reducing the 
dead zone) indicates that even this form of control is superior to 
the uncontrolled situation. An additional advantage is that as 
the degree of skid is reduced toward 0, i.e., as the relative 
velocity between the tire and road surface is reduced, the coef- 
ficient of friction becomes somewhat greater than that present 
when a more significant relative velocity is attained. This means 
that a greater value of deceleration is possible when the skid 
becomes extremely small. This feature of the coefficient of 
friction was not incorporated in these analog computer simulations 
due to the greatly increased complexity of the simulation that would 
be necessitated. 

Even though this control could theoretically provide good perform- 
ance if the dead zone were reduced sufficiently, there is a serious 
physical limitation that prevents this from ever really being possible. 
A little reflection will reveal that the lower the dead zone setting, 
the more rapidly will the corrective force be required to go on and 
off. As the dead zone approaches zero, the shape of the corrective 
force will appear as a high speed oscillation, and the practical lower 
limit for the dead zone will be determined by both the speed of response 
of the hydraulic or pneumatic lines and the physical strength of the 
brake system. 

Thus, from a practical point of view, this system appears to offer 
the advantage of a controlled skid, plus a shorter stop time than that 
provided by the double ramp control. The disadvantage is that a skid 
is almost always present, and the level to which it could be reduced 


would be severely limited by the physical response of the brake lines. 


50 


Neither control tried thus far appears to be superior to the 
other. This is verified somewhat by determining the performance 
factor for the bang-bang control with a dead zone. The performance 
factor (computed by taking the reciprocal of the product of the 
integral of the error and the difference between the stop time and 
the optimum time to stop) is determined to be 0.046. This is 
neither significantly better nor worse than any of the performance 
factors obtained for the double ramp control. 

The bang-bang control with a time delay was then simulated. 

The results for this simulation are presented in Fig. 23. A time 
delay of 0.1 sec. (problem time) was chosen for this simulation. 

The stop time is seen to be 5.93 sec., which is considerably worse 
than the stop time for the previously described bang-bang control. 
However, in this case the skid is eliminated for a portion of the 
total stop cycle, and the skid time is 3.92 sec. Superficially this 
appears to be an improvement over the previous bang-bang control, 
but closer examination reveals that very little, if any, betterment 
of the performance is gained by the incorporation of a time delay. 

Reference to Fig. 24, which is an expansion of a portion of 
Fig. 23, and was obtained by manual solution of simplified equations 
of motion, reveals that for a significant period of time during each 
error correcting cycle the velocity remains a constant. For example, 
from point 2) to (3) , where no skid is present, the velocity is 
also a constant. The time from (2) to (3) is 0.1 sec., and corresponds 
to the time delay built into the corrective force. During this time, 
the net force F on the brake shoe is zero, and thus the velocity 


remains constant. If the time delay is decreased, this control 
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approaches the performance of the bang-bang control with a zero 
dead zone. 

The bang-bang control with a time delay thus appears superior 
to a completely uncorrected braking situation, but appears some- 
what inferior to either of the two controls previously described. 

A performance factor for this control (for the specific set of 
parameters used for this simulation) is 0.025, which is sufficiently 
lower than the trend of the values obtained for the previous two 


controls to warrant the above conclusion, 
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7. THE TRIPLE RAMP CONTROL. 


The triple ramp control of Fig. 12 was the final type of 
corrective scheme to be implemented. A total of 100 separate 
computer runs were made, using a variety of values for the 
parameters K,> Ki K, and C.R. (correction ratio). Representative 
samples of the behavior to be expected from this type of corrective 
scheme are shown in Fig. 25 and Fig. 26. 

The most striking features about both Fig. 25 and Fig. 26 
are the low stop time, skid time, and magnitude of error present 
moc chic ftOmeoL G@ntrol. For example, Fig. 25, with Ky = 60000 
lbs/sec., K, = 40000 lbs/sec., K. = 1000 lbs/sec. and C.R. = 2.05 
has a stop time of 4.22 sec., skid time of 0.98 sec. and roughly 
1.9 ft/sec. as the degree of skid. This gives a performance factor 
of 0.493 which is higher than the maximum value attained for any 
control tried previously. This value of performance factor is of 
course affected by the values of the various parameters Ky Ko Kae 
and C.R., and does not represent the "best" value of P.F. that was 
obtained on the trial runs for the triple ramp control. Indeed, 
the term "best" is misleading here, since a setting of parameters 
resulting in the highest possible P.F. would sacrifice an adapta- 
bility to changing road conditions. This feature will be discussed 
more thoroughly later in this paper. 

Nevertheless, it can be concluded from the above factors and 
an inspection of Fig. 25 and Fig. 26 (along with the Author's 
analysis of the other 98 analog computer runs) that the triple ramp 
control provides more nearly optimal braking conditions than any 


previously described control. 
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Fig. 26, Triple ramp control response, Kj = 25000 lbs/sec., 


K, = 25000 lbs/sec., K, = 500 lbs/sec., C.R. = 2.05. 
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A closer examination of the actual behavior of this triple 
ramp control is appropriate at this point. The results of a 
manual solution of the equations of motion (with the simplifying 
assumptions mentioned previously for similar hand calculations) are 
presented graphically in Fig. 27. The parameter values assumed 
were K, = 50000 lbs/sec., Ky, = 80000 lbs/sec., K, = 1000 lbs/sec., 
and C.R. = 1.5. Also, the coefficient of friction w between the 
brake drum and brake shoe was assumed to be constant at 0.8 for 
the hand calculation (this was the only simplifying assumption). 
All other conditions were assumed the same as previously listed 
for other simulations. 

An examination of Fig. 27 indicates that from time = O sec. 
to the point marked oS , the behavior of the control and the 
vehicle response are precisely analogous to the double ramp form 
of control. This can best be seen by a comparison of Fig. 27 
with Fig. 17 over the region from time = O sec. to the point (5) : 
The same defect present in the double ramp control is also present 
in the triple ramp control, ive., from point (2) to point GB) 
the velocity is decreasing at a rate slower than the rate of 
deceleration when 2 skid is present. Also, from point (5) to point 
(6) this phenomenon occurs. However, at point (7) , when the down- 
ward slope of the corrective force B. abruptly changes from - Ky, 
to - K. the situation changes completely. While the slope of 
BF. =- Ko the rate of decrease of the velocity not only is signif- 
icantly greater than that possible during a skid, but also no error 


is present until the value of the deceleration again exceeds - 24.0 


Z : 2 , 
ft/sec. at point ; t this time the process repeats itself. 
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It is during this last phase of the skid correcting cycle 
(from point ) EO pOlme (8) in Fig. 27) that the inherent 
superiority of the triple ramp control is manifested. During 
this portion of the cycle the deceleration is gradually increasing 
as PF. slowly decreases. When the deceleration is too great for 
the prevailing road conditions a skid once again occurs, but prior 
to this occurrence of the skid the deceleration is more and more 
closely approaching the optimum deceleration. Additionally, if 
the road conditions are changing, for example, if the vehicle 
went from a patch of slick pavement to dry pavement, and the coef- 
ficient of friction between the tire and road surface were increased, 
then this form of control permits the deceleration to increase to 
the new optimum value. In other words, the triple ramp control 
is continually sampling the road conditions and permitting the 
braking force to adapt to changing road conditions. 

The desirability of this adaptability feature is evident. It 
is also apparent that for the invariant road conditions simulated 
for this problem, the most nearly ideal braking situation would 
have been achieved by setting C.R. to a value just slightly greater 
than one, and causing K. to equal zero. This would have permitted 
only two skids of small magnitude to occur, and for the remainder 
of the path the vehicle would have slowed down at a rate very nearly 
equal to the optimum deceleration. This would provide a very large 
value for the performance factor, which would of course be a misleading 
indication in this case, since the control would not respond to improved 
road conditions. It would, however, still respond to worsening road 


conditions, so there is even some merit in this proposal, 
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Referring back to Fig. 25, we note that the typical triple 
ramp shape of F. does not become evident until after the occurrence 
of the fourth error signal, whereas in Fig. 26 the triple ramp 
shape appears as expected after the occurrence of the second error 
signal. Recall that the form of the coefficient of friction wu 


between the brake drum and brake shoe is given by 


tie= 158 5= 5. 00e@ - 60) for w > 60 rad/sec. 
Reference to the W versus time plot of Fig. 25 reveals that both 
the first and second errors occur while the rotational velocity is 
above 60 rad/sec. Thus in this region the value for » is constantly 
changing in accordance with the above relation. Now, when the first 
error signal (or skid) is sensed the corrective force BF. immediately 
increases until the skid is eliminated, and then decreases until 
the skid occurs once more; this value of F. at which the skid occurs 
is stored by a storage device. Upon the occurrence of the second 
skid, PF. increases again until the skid is once more eliminated. 
Ideally the value of F. should drop to some multiple of the stored 
previous value of F. (the multiple being established by C.R.). How- 
ever, the velocity of the vehicle is decreasing continually during 
this process and w is thus increasing. This means that after the 
occurrence of the second skid, the net torque on the wheel for a 
given value of F. (e.g., for a value equal to C.R. times the previously 
stored value of F.) is now greater than it would be if the value for U 
were not increasing. Thus a2 third skid occurs before F. can drop to 
the value at which the break point occurs. However, after the occur- 
rence of the fourth skid, w% is constant, and thus the typical triple 


ramp effect is observed. 
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In Fig. 26 the above phenomenon does not occur. Note that 


in this figure the value of K, is only 25000 libs/sec. whereas in 


3 
Fig. 25 K, had a value of 60000 lbs/sec. This means that the 


3 

error will be corrected more slowly, and thus the velocity will 
have decreased to a level where w™ is a constant by the time the 
second error has been eliminated. Thus the coefficient of friction 
will not have changed sufficiently between the time when the value 
of F. was stored and the time when B. reaches its break point to 
cause another error signal to occur. Note, however, that if C.R. 
had been set sufficiently low (e.g., slightly above 1.0) then the 
Same phenomenon that occurred in Fig. 25 would also have occurred 
in Fig. 26. Conversely, had C.R. been set sufficiently high, this 
phenomenon would not have occurred in Fig. 25. 

Figure 28 presents a plot of stop time versus Ko with @gRe 
as a parameter. The data points on this and subsequent plots 
were obtained from the graphical output resulting from the various 
analog computer runs made while simulating the triple ramp control. 
Consequently, the accuracy of the individual data points is only 
of the order of a few percent. However, this degree of accuracy is 
sufficient to demonstrate the overall performance characteristics 
of the triple ramp control quite satisfactorily. From Fig. 28 we 
can see that as K. increases the stop time increases also. However, 
the rate of increase, i.e., the slope of the curve, is observed to 
be a variable, even though the overall trend is upward. This over- 
all upward trend is expected in view of previous comments concerning 
the adaptability feature inherent in the triple ramp control. A 


decrease in stop time must be achieved at the cost of a decrease in 
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Fig. 28, Stop time versus K. with C.R. as a paraneter 


for the triple ramp control, 
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the sensitivity of the control to improving road conditions, all 
other parameters being assumed constant. Thus, as K. increases, 
the vehicle goes into a skid mode more frequently, and the overall 
trend of stop time with K. would be an upward one. 

However, to satisfactorily account for the leveling off of 
the stop time versus KR, chart at localized regions, we must investi- 
gate the detailed behavior more closely. Figure 29 is a plot of 
both stop time and skid time versus K. for one specific set of para- 
meters. The numbers adjacent to the skid time curve represent the 
total number of skids present in the particular graphical output 
from which each specific data point was obtained. For example, 8 
skids were present in the analog computer output obtained with a 
setting of K, = 60000 lbs/sec., Ky, = 25000 Ibs/sec. ,9G.R. = 2205. 
and K, = 1500 lbs/sec. For K, = 2500 Ibs/sec. and K, = 3000 lbs/sec., 
10 skids were present and the skid time was thus the same for each 
of these values. However, Fig. 29 also reveals that the stop times 
for each of these values of K. were the same. 

In a situation where the same number of skids is observed from 
one value of K. to another, then the skid times will obviously be 
identical, assuming that K, (which directly determines the magnitude 
of each skid and the length of time it is present) is the same for 
each trial. The only factor contributing to a difference in over- 
all stop time is thus the difference in time between adjacent pairs 
of skids (e.g., the time from point (5) to point stl rei) 

For a larger value of Ko the time between pairs of skids is decreased. 
Now, for example, when 10 skids occur for a value of K_.= 3000 lbs/sec. 
the time between adjacent pairs of skids is shorter than the time 


between the adjacent pairs of 10 total skids obtained for K. = 2500 
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Fig. 29. Skid time and stop time versus K, for the triple 
ramp control, K3 = 60000 lbs/sec., K, = 25000 lbs/sec. 


C.R. = 2,05. 
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lbs/sec. Thus the final pair of skids occurs earlier when K. - 
3000 lbs/sec., but due to a larger value of K. a larger value of 
deceleration is attained after the tenth skid and prior to the 
vehicle coming to a complete stop. There is also a longer period 
of time after the tenth skid occurs, during which the slope of 

PF. = - K. for the case where K, = 3000 lbs/sec. These two 
factors serve to account for the leveling off of the stop time 
versus K. curves observed for Fig. 28 and Fig. 29. It is also 
possible that a slight local decrease in total stop time might 
occur for an increase in KR, under the proper combination of 
control parameters and road conditions. This is observed in Fig. 
28 in one region of the curve obtained for C.R. = 2.74. 

Despite the above observations, the conclusion to be drawn 
from Fig. 28 and Fig. 29 is that the overall effect of an increase 
in K. is to increase the stop time of the vehicle, and to degrade 
the overall performance of the corrective system. A compromise 
must thus be reached between adaptability to improving road 
conditions and decreasing the total stop time. Designing for a 
worst case situation, a satisfactory value for K. would thus be 
K. = 500 lbs/sec., as this would permit reasonably rapid response 
to improving road conditions while still providing very good braking 
performance under nearly static road conditions. 

The next parameter to be investigated is the correction ratio 
(C.R.). Figure 30 is a plot of the stop time versus the correction 
ratio, with K. as a parameter. The overall trend as determined from 
Fig. 30 is that the stop time tends to decrease with an increase in 


the value of C.R. Each curve in the family of Fig. 30 demonstrates 
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notable local exceptions to this trend, however, and, as will 
be pointed out, excessively large values of C.R. will materially 
lengthen the stop time. 

The overall downward trend of the stop time versus C.R. curve 
is best understood by referring to the inset of Fig. 27, where a 
pictorial description of the triple ramp control is presented. 

As the value of C.R. is raised this raises the break point at 
which the transition of the slope of B. from - Ky toe K. occurs. 
This has the effect of lengthening the time between the occurrence 
of adjacent skid pairs, since the time that B will take to drop 
to the level where a skid will occur once more increases in direct 
proportion to the height of the break point. It is thus seen that 
increasing C.R. has the same net effect as decreasing the value 

of Ka» with certain minor restrictions. 

It would be expected that, as C.R. is increased to obtain a 
set of braking histories, those output traces which have the same 
number of skids present for differing values of C.R. would behave 
in a manner Similar to that observed when K. was the variable 
quantity (e.g., in Fig. 28 and Fig. 29). Referring in Fie, 30 to 
the curve obtained for K. = 500 Ibs/sec., we see that the stop time 
increases with increasing C.R. when the total number of skids present 
in each braking history remains a constant four (this is also observed 
for K, = 1000 lbs/sec., for six skids observed). The reason for 
this is very simple. As C.R. is raised this "pushes" the adjacent 
pairs of skids farther apart and thus for a larger portion of the 
total stop time the vehicle is decelerating at a lower rate than if 


C.R. were set lower. Thus if the total of Skids Se wie ren oimed tic 


68 


same, the control which provided the greater deceleration over 
the entire skid free path would also provide a shorter stop time. 
This is the control with the lower C.R. 

Note that C.R. is also closely linked to the adaptability of 
the system to improving road conditions, since for a large value 
of C.R. the time required for Bo to drop (at a rate of - K.) to 
a level where the road is "tested'' for better conditions is 
lengthened. Additionally, C.R. cannot be increased to a very 
large value without greatly extending the stop time. The minimum 
number of skids possible (starting with any significant initial 
velocity) is two. Once C.R. has been increased to a level where 
only two skids are present, then any further increases in C.R. will 
cause a rapid increase in the stop time for the reasons listed in 
the preceding paragraph. However, long before the skids were 
reduced to two in number, the stop time would also have increased 
due to the low value of deceleration that would be present over an 
excessively large portion of the vehicle's trajectory. 

In retrospect, it appears that despite the downward trend of 
stop time with increasing C.R. presented in Fig. 30, any advantages 
to be gained by very large values of C.R. are outweighed by the 
disadvantages. On the basis of information presented here, a "safe" 
value for C.R. would appear to be about 2.0. 

Figure 31 and Fig. 32 present plots of both stop time and skid 
time versus C.R. for two different sets of parameters. They are 
included primarily for information, the most significant factor 
about either one being that where the curve representing skid time 


becomes flat (implying a constant number of skids were observed for 
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succeeding braking history plots) the curve for stop time acquires 
a positive slope. This serves to reinforce the arguments previously 
advanced, 

The behavior of the braking performance as Ry is varied is 
also of some interest. Figure 33 illustrates the behavior of stop 
time and skid time as K, is increased. Notice that the tendency 


5 


exhibited is for the stop time to decrease with increasing K This 


3° 
is the expected result, since the larger the value of 3 the more 
rapidly is any error or skid corrected back to zero, and thus a 
larger portion of the braking path will be in a controlled mode of 
skid free braking. 

Referring back to Fig. 27 once more, we note that the time 
from the commencement of braking action (t = 0") to point , and 
also from point (4) to point (5) , 1s dependent solely upon the value 
of KR, whereas the time from point Q to point (4) and from point (S) 
to point depends upon Ki, and Nat Thus as K, is increased the 
total skid time and magnitude of each individual skid are reduced, 
This has the effect of causing each successive skid pair to occur 
earlier in time, even though the time between skids is unchanged. 
As long as the total number of skids present is unchanged as K, is 
increased, an improvement will always be observed in the stop time. 
However, in the event that one value of K., causes a certain number 


3 


of skids, and a higher value of K, is just sufficient to increase 


3 
the number of skids, then a slight local increase will be observed 
in the stop time. This phenomenon is seen to occur in Fig. 33. 
The explanation for this is as follows. When K, = 50000 ibs/sec. 


the time remaining between the correction of the sixth skid and the 


complete stopping of the vehicle is insufficient for the deceleration 
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Fig. 33. Stop time and skid time versus Ky for the triple 
ramp control, K, = 25000 lbs/sec., K, = 1000 lbs/sec. 
C.R. = 2.05 
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to have increased to that level necessary to cause another skid. 
However, for K, = 60000 lbs/sec. the seventh skid occurs just 
prior to the vehicle coming to a rest. During this seventh skid 
the vehicle deceleration is only - 16 fe) scene rather than the 

- 24 reeeer that was experienced just prior to the skid. It is 
this lower value of deceleration present just prior to stopping 
that causes the stop time to increase as K, is increased from 
50000 lbs/sec. to 60000 lbs/sec. 

The above phenomenon is just restricted to the transition 
points where the number of complete skids increases from one value 
to another. However, the decreased skid times and magnitude of 
skids caused by an increased value of KR, indicates that the greatest 


overall advantage is to be gained by making K, as large as possible, 


3 
regardless of any localized increases in stop time caused by this 
POUNCY . 

Finally, the system behavior with changing values of Ky is 
investigated. Figure 34 indicates that the stop time decreases as 
the value of Ky is increased, but that the skid time goes up with 
increasing Kye Reference to Fig. 27 again facilitates understanding 
this behavior. From Fig. 27 it can be seen that the time lapse from 
point (2) to point (4) and also from point GS) to point @D is depend- 
ent solely upon the value of Kye Large values of Ky, will cause the 
time intervals between the points to become small, while leaving 
time intervals along other portions of the trajectory unchanged. 
Since it is during the time interval from point (5) to point @, 


for example, that the deceleration of the vehicle is even less than 


that present while in a skid, any factor which will shorten the time 
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from ©) to @ will also shorten the time the vehicle spends in 
the undesirable region from (5) to 6) » This will naturally 
cause a shorter stop time, since the vehicle will be decelerating 
at a high rate over a larger proportion of its stopping path. 

The increase of skid time noted in Fig. 34 is once again 
related closely to the total number of skids present. As the time 
' from (2) to (4) or from (6) to @ (referring once again to Fig. 27) 
is reduced, this tends to shorten the time between the occurrence 
of adjacent skids. If the total number of skids remains the same 


while K, is being increased, the skid time naturally remains the 


os 
same (to fully appreciate this, one must recall that the duration 
and magnitude of each individual skid is controlled by K, only, 
which is not being varied here). However, as the skid occurrences 
become closer together in time, a point is reached where the total 
number of skids encountered in the stopping history increases, and 
thus the skid time will also increase. This increase in skid time 
would not be a serious consequence, however, if Ky were sufficiently 
high to make the duration and magnitude of the individual skids very 
small. 

Thus it might be concluded that the larger the value of K,, 
the better the performance of the entire system, since the decrease 
in stop time would more than offset any increase in skid time, partic- 


ularly if K, were also set very high. Caution must be exercised here, 


3 
however, because it is the nature of the triple ramp control that 
the value of B at which every second error occurs must be marked 


and stored. This means that if Ky, is too large, then difficulties 


would be encountered in marking the point at which each second error 
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occurs, and additionally in reversing the trend of Fe (i.e., in 
changing from a slope of - K, to ate Ky )- Obviously a slope of 
- infinistyestor - Ky, would be totally unsatisfactory. Thus, a 
very large value of Ky, would indeed provide the best performance, 
but this value must not be so great as to create excessive dif- 
ficulties in measuring and storing the point at which each second 
error occurs. 

In summary, the following qualitative comments concerning 
the choice of the values for the various parameters of the triple 
ramp control are appropriate. K, should, with no qualifications, 
be made as large as possible. Ky should also be as large as possible, 
short of being a negative step function, and subject to the restric- 
tions mentioned in the preceding paragraph. K. should be relatively 
low assuming that C.R. is not too large. A value for C.R. of about 


2.0, and for K_ of 500 lbs/sec. was determined to be satisfactory 


for the particular system simulated here. 
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8. PROPOSED IMPLEMENTATION OF THE TRIPLE RAMP CONTROL. 


In view of the obvious superiority of the triple ramp form 
of control over any of the other controls examined, a proposal 
for the physical implementation of this control is advanced here. 
It must be assumed that a linear accelerometer of high accuracy 
and linearity over the range from 0 pee to 32 ‘ual SEs is 
available. A high gain, broad bandwidth operational amplifier is 
necessary to provide the integrator used in conjunction with the 
linear accelerometer (see Fig. 1 for the block diagram of the 
system). Additionally, an accurate tachometer to measure the 
rotational velocity of the wheels is assumed to be available. 
Thus, all elements necessary to accurately register the occurrence 
of a skid are assumed to be available. Our attention will be 
directed toward the actual mechanism by which the triple ramp 
control is generated, and by which the pressure applied to the 
brake cylinder is modified. 

Figure 35 illustrates the means by which it is proposed to 
implement the triple ramp control. An explanation of Fig. 35 
follows. 

A force which is applied to the brake pedal is transmitted to 
the brake drum via the hydraulic lines. Until a skidding situation 
occurs Vo is held equal to Nye by means of the ramp generator and 
limiter. When VW equals Vee full braking action is applied and 
no corrective action is occurring. 

The first time a skid occurs, e (error) becomes greater than 


zero, and switches (1) and @ throw in the indicated manner. 
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Proposed implementation of the triple ramp control. 


Capacitor C, discharges rapidly through resistor Ry (a very small 


3 
resistance). Us thus drops very rapidly, permitting the braking 
force to be decreased toward zero; this step is an exponential 
decay with a very small time constant, but it approximates the 
first upward ramp of B. (shown with a slope of Ky). 
When the error is reduced to zero (skid eliminated) switches 
@) and (2) throw once more. The voltage on C, then starts to climb 
with a slope of K,, thus causing the braking force to increase once 
more. During this first part of the cycle a portion of as is being 
monitored by C, (thru the tap on Rj). When the braking force has 
increased sufficiently to cause skidding to occur once again, 
switches @ and 2) throw once more and Nie drops rapidly, decreasing 
the brake force as before, Additionally, however, switches G3) and 
(4) are now thrown, and capacitor Cc. has a charge proportional to 
the force required to cause skidding (say .9 x ve to cause a skid). 
When We falls sufficiently to eliminate the skid, e = O again, 


and the ramp generator charges C. upward once again. When Vo is 


> 

just equal to the charge on C, CG, ox Ve to cause a skid, for example) 

V becomes greater than zero, and switch 6) throws as indicated. 
Eese 


This puts a Large capacitance (C,) in parallel with C, and causes 
the ramp slope to change to a much shallower slope (i.e., to change 
from slope = Ky to slope = Ko). ve is thus caused to approach the 
value that creates a skid much more slowly. 

When (or if) a skid occurs once more, the entire process 
described above repeats itself. 


A sketch of the force acting on the brake cylinder is shown in 


Bao 


80 


TIME CONSTANT =RaCz 





Fig. 36. Force acting on the brake cylinder during occurrence 
of one skid pair, if using Proposed Implementation. 

Note also that the sensing of the first, third, etc., skids 
or second, fourth, etc., skids can be accomplished by an R - S 
flipflop connected to the output of a comparator which would be 
used to indicate the presence or absence of a skid. This was 
done on the analog computer simulation with very satisfactory 
results. 

Figure 37 is a blow up of the pressure controller shown in 
the inset of Fig. 35. This is a standard form of control valve, 
and is not original with the Author.[1] A brief explanation of 
the workings of a standard power brake system, for which this 
particular form of implementation is ideally suited, is appropriate 
here. 

In a power brake unit, when the brakes are applied, a valve 
to the atmosphere is opened on one side of a piston, and the other 
side is connected to a vacuum. The atmospheric pressure forces 
the piston closed, and causes pressure to be transmitted to the 


brake cylinder via a hydraulic line. 
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Fig. 37. Detail of pressure controller suitable for use 
with a power brake system. 

The above controller will cause a disruption of the pressures 
on the power cylinder if the valve is rotated from the position 
shown. In the position shown, the vacuum line is blocked off thus 
poe cine full atmospheric pressure to the piston; any other 
position of the valve would lessen the pressure on the piston, 


thus reducing the braking action. 
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9. CONCLUSIONS. 


Of the three major types of controls simulated, the triple 
ramp control proved to be the only one lacking major defects, in 
that it reduced the degree of skid to a very low level while still 
providing a very short stop time. Nevertheless, the double ramp 
control was determined to provide a definite improvement over the 
completely uncontrolled situation, as did the bang-bang controls. 
However, for the double ramp control a slightly increased stopping 
time was the penalty paid for obtaining a lower degree of skid 
while stopping the vehicle. The bang-bang control with a dead 
zone proved unsuitable due to the presence of a slight degree of 
skid during the entire braking path of the vehicle. Additionally, 
the level to which the skid could be reduced would be severely 
limited by the response characteristics of the brake lines. The 
bang-bang control with a time delay, while still superior to the 
uncorrected situation, appeared to be considerably inferior to any 
other form of control simulated. 

The triple ramp control, however, appears to provide very 
nearly optimum braking conditions in that the degree of skid and 
the skid time are both minimized, while still providing a considerably 
shorter stop time than would be obtained from any other form of 
control. A trade off between stop time and adaptability to improving 
road conditions is necessary though, and the proper balance between 
these two factors would be determined by both the type of vehicle 
involved and the anticipated road condition. However, greater 
emphasis should probably be placed on minimizing skid time and stop 


time by having a relatively low value for Kae 
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The problem of possible oscillation introduced by the triple 
ramp control should never arise, since it is essentially a negative 
feedback system and should be inherently very seen ete ly 
unusual road conditions would probably be necessary for any oscil- 
lations to occur; however, an investigation of this subject has 
not been undertaken. This would possibly be a fruitful area for 
further study. 

The response time required for the mechanical components of 
the proposed implementation of the triple ramp control (Fig. 35 and 
Fig. 37) would certainly be easily attainable. This may be seen 
from Fig. 25 (which is a "typical" response curve), where examina- 
tion of the F. versus time curve indicates that the pressure control- 
ler would go through the high speed cycle at the rate of only about 
4 c.p.s. This serves to give an order of magnitude approximation 
to the response characteristic necessary. 

Throughout this entire thesis, Pea has been assumed to be a 
constant, implying a panic stop situation with maximum force applied 
to the brake pedal. If, during the course of slowing the vehicle, 
the operator decreased the pressure applied to the brake pedal, no 
effect would be noticed in the operation of the automatic control, 
as long as the value of Fa were still sufficiently large to cause 
a skid (if acting by itself). If, however, Ee were decreased below 
the value needed to cause a skid (if acting alone), then the control 
would permit the net force F to rise to the new value of Fea (either 
at the rate of Ky or Kae depending upon which part of the cycle was 
acting at the instant of decreasing the value of Few) and this value 


of F would be retained until a skid once again occurred. If the 
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applied brake force Pan were increased, the action of the triple 
ramp control would be unaffected, as the control would still 
continue to seek the optimum net force F. 

A matter of possibly more concern is the characteristic of 
the linear accelerometer and of the integrator used to indicate 
the existence of a skid. The proper operation of this skid detect- 
ing scheme is basic to the operation of any skid correcting device 
proposed in this thesis. The mounting of the linear accelerometer 
would be very critical, and a position on the axle housing, in 
order to minimize ''nose dive'' effects, would be the most nearly 
ideal location. However, no serious theoretical limitations on 
the effectiveness of this form of skid detecting scheme present 
themselves, and this would appear to be the only practical means 
of detecting a skid prior to actual wheel lock up. 

A final comment concerning the analog computer on which the 
brake system and all controls were simulated is appropriate at 
this time. The COMCOR Ci 5000 Analog Computer is in actuality a 
form of hybrid computer, in that extensive logic circuitry is avail- 
able for use in conjunction with standard integrators and operational 
amplifiers. Simulation of this type of problem would be difficult 
1f not impossible on a more conventional analog computer lacking 
the logic circuitry, high speed switches and overall versatility 
of this particular computer. A simulation of the type described 
would be possible on a digital computer utilizing one of the 
languages available for this purpose, such as DSL-360, etc. How- 
ever, in order to obtain the extensive data needed to determine the 


performance characteristics of the various controls, large quantities 
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of very expensive digital computer time would be required. On 
the analog computer, however, almost two hundred different runs 
were made consecutively in a relatively short (and inexpensive) 
time period by merely varying potentiometer settings between runs. 
In this case, the accuracy of a digital computer was unnecessary, 
and the (hybrid) analog computer proved to be an exceedingly 
valuable design aid. This thesis thus illustrates graphically 

the continuing need for increasingly more sophisticated analog 


computers, such as the COMCOR Ci 5000. 
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APPENDIX A 


Determination of the maximum possible deceleration of a vehicle. 
The maximum deceleration possible for a vehicle is determined 

by the coefficient of friction between the tires and the road 

surface. Figure Al represents an idealized vehicle traveling 


in the plus direction with a velocity X, and a deceleration X, as 


indicated, 









Weight = mq x , 
x ——> 
ONG 
—<— Fr 


Fig. Al. Idealized vehicle. 


In Figure Al, Fe represents the force of friction acting on 


the vehicle, and is the only external force causing the vehicle 
to decelerate. The weight, mg, is the normal force acting on 


the ground, so the value of Fe is given by 


Fe= pmg , (Al) 


where WY is the coefficient of friction between the wheel and the 
road surface. The coefficient of friction is a variable, which 
has its greatest value (bo) just before a skid occurs, and then 
drops abruptly to a lower value Ly when a skid begins. 


Thus the maximum value of Fe is 


eG) : ens : (A2) 
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Note that corresponds to the coefficient of kinetic friction, 


RK 


whereas Le corresponds to static friction. 


Summing the forces in the X direction, we obtain 


mX = - Fe : (A3) 
Substituting Eq. (A2) in Eq. (A3) and solving for ey? we 
obtain for the maximum possible deceleration 
- eae “E68 - ee 
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tien of C 


APPENDIX B 


ls used 


The analog computer elements used in the simulations described 


SYMBOL 
"i ! 
Xp | X out 
< 10 
I.C.=Xo 


x 


XK 
fu a 
: 
a 
E 
ct 


xX 

Xrest > K 

Xout 
.___J 
Xtrest >K 

K, — Xout 
XK C 
Xe C 


in this thesis are indicated below. 


DESCRIPTION 


Sign changing summer amplifier. 


Xout = — (XK, + Xo + 10 xs) 


Sign changing integrator. 


Gs. 


Xeue= —| f (xX, +1OXa) dt + Xo 


To 


Class III multiplier. 


Xout = X, Xo 


Potentioneter. 


Xout=K,X; » K,< | 


Double pole, double throw relay. 


Xone = X, ) For Xrest << K 


Xout = Xa for Xrest > K 


High speed digital/analog switch. 


Xout = O yfor Xrest<k 
Xout = Xi) for Xrestr> K 


Comparator. 


C = logical 1 and C = logical 0, 
when xX; + x5> 0. 

C = logical 0 and C = logical 1, 
when x; + x5<0. 

Comparators are used to throw both 


the DPDT relays and the D/A switches. 
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26. 


——.. 


Xi 


resis 





TESTs> Ko 


Xout 


Tie point. 


Used in conjunction with two or more 
D/A switches which are thrown so that 


only one switch is ON at a tine, 
Koce a x Or Xe 


Switching accomplished by NAND logic, 


utilising only one D/A switch, 


Xout a A J when 
CREO | -aaweond “Ch EGiks> Ke). 
Also , Xout =O, When 


CMe Koon (TE Siem Ka) - 


ou 
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